The leading Ruelle resonances of chaotic maps 
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The leading Ruelle resonances of typical chaotic maps, the perturbed cat map and the standard 
map, are calculated by variation. It is found that, excluding the resonance associated with the 
invariant density, the next subleading resonances are, approximately, the roots of the equation 
z 4 = 7, where 7 is a positive number which characterizes the amount of stochasticity of the map. 
The results are verified by numerical computations, and the implications to the form factor of the 
corresponding quantum maps are discussed. 
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The statistical charectaristics of the quantum mechan- 
ical spectrum of a system, which has a classical analog, 
are related to its underlying classical dynamics. For ex- 
ample, in chaotic systems, this relation is revealed by 
Gutzwiller's trace formulafil which expresses the density of 
states in terms of a sum over the classical periodic orbits 
of the system. In disordered systems, spectral properties 
of the quantum system are expressed in terms of the dif- 
fusons and the Cooperons diagrams associated with the 
classical diffusive modes of the systema. 

There is, however, a difference between the periodic 
orbit picture and the diagrammatical approach of disor- 
dered systems. The first uses the individual trajectories 
of the systems as the basic ingredients for the semiclas- 
sical expansions, whereas the second approach employs 
spectral properties of the evolution operator of classical 
distribution functions. The latter approach is, therefore, 
more suitable for a field theoretic treatment. 

Recently an extension of the field theoretic formal- 
ism, from disordered systems to general chaotic ones, has 
been proposed!! The analog of classical diffusion, in the 
generalized case, is the coarse grained Liouville dynam- 
ics of distribution functionscl in the limit of zero coarse 
grainingfl. The corresponding evolution operator is the 
Frobenius-Perron (FPX operator, and its eigenvalues are 
the Ruelle resonances!! These eigenvalues describe the 
decay of classical correlation functions of smooth observ- 
ablcs. 

The association of spectral properties of quantum 
chaotic systems with their classical spectra calls for a 
study of the Ruelle resonances of general chaotic systems. 
The central question is whether chaotic systems can be 
classified into equivalence classes according to their clas- 
sical spectrum? An obvious example of such a class is 
systems where the effective dynamics is diffusive. For 
instance, the kicked rotor map exhibits diffusion in mo- 
mentum spaed!, and the stadium billiard, in the limit of 
small distance between the semicircles, is diffusive in the 
angular momentum spacen 

Are there other universality classes? In principle, there 
should be. Diffusion is not the only type of universal dy- 
namics of chaotic systems where correlations are rapidly 
lost. Other possibilities can be, for instance, various 



types of Levy flightsEI. 

Yet, the lack of approximation schemes for calculat- 
ing Ruelle resonances, obstructs the characterization of 
the effective dynamics of general chaotic systems. The 
physical approach for calculating Ruelle resonances is to 
project the dynamics, from the full phase space, down 
to the manifold on which the dynamics is slow, and to 
construct an equation for the probability density on this 
manifold. For example, in disordered systems, the fast 
and the slow manifolds are the momentum and the real 
space respectively. Projecting the dynamics down to real 
space yields the diffusion equation, whose eigenvalues 
constitute the leading Ruelle resonances of the system. 
In trying to apply a similar procedure for general chaotic 
systems, one encounters the problem of identifying the 
slow and the fast manifolds. These can be complicated 
functions in phase space. 

In this paper we choose a different approach. We will 
construct a simple variational method for calculating the 
leading resonances, and— apply it to two generic-maps: 
The perturbed cat mapEj and the standard mapQ. The 
results indicate that, indeed, the leading Ruelle reso- 
nances of these maps share some universal features (dif- 
ferent from diffusion) . Namely, the configuration of the 
leading resonances in the complex plain is similar. Ex- 
cluding the eigenvalue associated with the invariant den- 
sity, Zo = 1, the next subleading resonances are, approx- 
imately, the roots of the equation: z 4 = 7, where 7 is a 
real positive number, smaller than unity, which decreases 
as the map becomes more stochastic. It will be shown 
that this configuration of resonances is significant for the 
spectral statistics of the corresponding quantum maps. 
For example, in the form factor, it leads to a suppression 
of the nonuniyersal corrections to the results of random 
matrix theoryB (RMT). 

In introducing our ideas it will be instructive to con- 
sider a seecific example. We begin with the perturbed 
cat maplMj which is an area preserving map from the unit 
torus to itself: 



•£71+1 2x n 4- y n 

y n+ i = 3x n + 2y n + f(x n+1 ) 



(Mod 1). 
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Here (x n ,y n ) are the phase space coordinates at discrete 
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graining s — 0.001. The dark and the bright regions corre- 
spond to the minima and maxima of respectively. This 
function follows the unstable manifold of the map which is 
shown on the right panel. 



time n, and 

k 

f(x) = — [cos(27rx) — cos(47ra;)] 

is a perturbation function. This function is smooth and 
periodic on the torus, therefore, the map is continuous, 
and diffejeentiable. At k — Or-the map, known as Arnold's 
cat mapt3, is highly chaoti<£3. Namely, time correlation 
functions decay faster than exponential. This behavior is 
related to the presence of hidden symmetries of number 
theoretical origiiJlJ. The role of f{x) is to break these 
symmetries. As a result, at finite values of k, relaxation 
is exponential. It has been also proved that for k < 0.11 
all the fixed points of (Q) are in one to one correspondence 
with the original cat map. Consequently, all orbits are 
unstable-and the system can be characterized as "hard 
chaotic"lil 

The classical propagator of the map, for one time step, 

is 

(x',y'\U\x,y) = 8 p [x'-2x-y]8 p [y'-3x-2y-f(x% (2) 

where 8 p (x) is a periodic 8- function on the torus. Due 
to Liouville's theorem, the above propagator is unitary 
when acting on Hilbert space of square integrable func- 
tions. Within this space the dynamics is reversible, how- 
ever, as soon as the dynamics is coarse grained, it be- 
comes irreversible and phase space densities relax to the 
invariant density. 

A simple way of coarse graining the dynamics is to give 
the (5-functions in (||) a finite width, for instance, by the 
replacement 

S p (x) -> V -^e-^-^, 

n v 

where s is the coarse graining parameter. The Ruelle 
resonances of the map are obtained by diagonalizing the 
propagator at a finite value of s, and taking s to zero 



after the diagonalization. It can be shown that this pro- 
cedure is equivalent to the introduction of noise into the 
equations of motion (|]), calculation of the eigenvalues av- 
eraged|Over the noise, and finally taking the limit of noise 
to zeroll3. The resulting resonances have absolute values 
smaller than unity, except for the resonance associates 
with the invariant density, Zq = 1. 

Since the coarse grained evolution operator is not uni- 
tary, the left eigenfunctions of U are not the complex con- 
jugates of the right eigenfunctions. Instead they can be 
viewed as the right eigenfunctions of the propagator as- 
sociated with the inverse map (i.e. the propagator taking 
the system one time step backward) . Thus the eigenvalue 
equations are: 

U\tf*) = z a \lg), (^\U=(^\z a , (3) 

where \4>^) and {ipa\ are * ne right and the left eigen- 
functions associated with the resonance z a . According 
to Ruelle's theoremQ, these eigenfunctions become distri- 
butions in the limit of zero coarse graining. The distribu- 
tions associated with the left and the right eigenfunctions 
follow the stable and the unstable manifolds respectively. 
Nevertheless, they satisfy the normalization condition 

In the left panel of Fig. 1 we show an example of the 
right eigenfunction, ip^, of the perturbed cat map for 
k = 0.5 and s — 0.001. As can be seen, this eigenfunc- 
tion follows the unstable manifold depicted on the right 
panel. As the coarse graining parameter, s, reduces, the 
eigenfunction becomes more singular, in accordance with 
Ruelle's theorem. 

In what follows we shall calculate the leading Ruelle 
resonances of the map ([!]) using variational approach 
which takes into account the structure of the eigenfunc- 
tions illustrated in Fig. 1. For this purpose we construct 
the functional 

F=(L\U\R)-z(L\R), (4) 

where \R) and (L\ are the right and the corresponding 
left trail functions. The leading Ruelle resonances are 
obtained by variation!!] of F with respect to (L\. 

The success of a variational calculation depends on the 
extent to which the variational wave functions capture 
the physics of the problem. In our case, one should con- 
struct left and right trail functions which follow the stable 
and the unstable manifolds respectively. We will use the 
dynamics itself to generate these functions by repeated 
applications of the map on a smooth initial state. Thus, 
our right trail function takes the form 

oo 

\R)= ]T'A„|n), with \n) =e 2 ^^y\ (5) 

n— — oo 

where A n are amplitudes to be determined by varia- 
tion, and the phases 4>n(x, y) are defined by the recursion 
equation: 
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<f>2(x, y) = ^ x {2x -y + f(x),2y-3x- 2f(x)) 

with the initial phase (f>f (x, y) = x. We define <j)^ n to be 
—<j>ni an d the prime ( ') indicates that the sum does not 
include the n = term. 

The above sequence of phases, where the first few of 
them are 
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6f 
b% = 2x 
bf = 7x 

6f - 



y + f(x), 

iy + Af(x) 



mi 



(6) 



= 26a; - Ihy + 15f(x) + 4/(^) + /(</#), 



is constructed by application of the map (Q) on the ini- 
tial smooth state, e t2 ^ x . Therefore, |2),|3),--- become 
increasingly smooth along the unstable manifold and sin- 
gular along the stable one. 

In the same manner we construct the left function by 
repeated application of the inverse map, namely 



(L\= 



(m\ with (m| = e - 2 ^™(^) ; (7) 



where the left phases are defined by the recursion relation 

<fa{x,y) = ( f>i_ 1 (2x + y,3x + 2y+ f(2x + y)) 

and (f>i (x, y) — x. 

There are two simplifying features characterizing \R) 
and (L\. First, the operation of the propagator, U, on 
these functions is simple since, by construction, 

U\n) = |n + sign(n)), and (m\U = (m + sign(m)|. (8) 

The second feature is that (L\ and \R) are orthogonal 
to a constant function. This way we eliminate the lead- 
ing Ruelle resonance, Zq = 1, from our problem, since 
the invariant distribution is constant in phase space. 
The proof for this orthogonality is straightforward. Let 
|0) = (0| = 1 denote the invariant distribution. Then 
(0|n) = <0| t^"— 1 1 1> = (0|1> = 0, where we have used the 
definition of the invariant distribution: (0\U = (01, and 



that II) = 



is orthogonal to a constant. Now, from 



(g) it immediately follows that (0\R) = 0. Similarly it is 
straightforward to prove that also (L|0) = 0. 

Variation of F, given by (||), (||) and (0), with respect 
to A* m yields the eigenvalue equation: 



Det [(m|n + sign(n)) — z(to|?t,)] = 0. 



(9) 



As it stands, this equation is not simpler to solve than 
the original problem. However, it can be considerably 
simplified if the matrices (m\n) and (m\n + sign(n)) can 
be truncated to a small size. As we show below this is 
indeed the situation. For this purpose, it is instructive to 
notice few properties of the matrix elements (m\n). The 
first one is that 



(m\n) = (m — sign(m)|7i + sign(n)), 



(10) 



where \m\ > 1. This is an immediate consequence of 
(m\n) = (m|?7 _1 [/ri), and the relations (|J). The second 
property of the matrix elements comes from the defini- 
tion of the phases with negative indices, 



{— to — n) = (to n) 



and 



(to — n) — (— to \n) 



(11) 



From ( |l(i| ) it follows that all matrix elements (±to| ± n) 
can be reduced to integrals of the form 



(±1 



±0 = / 



dxdye 



27ri[=Fx±0f (x,y)] 



where i = |m| + |n| — 1. Notice that the above integrals 
decrease rapidly with I. This follows from the nature of 
the phases cf>f (x,y). Namely, as I increases, the phases 
acquire a stronger dependence on x and y, see e.g. (](]). 
Thus (m\n) is small for large \n\ + m\, and the signifi- 
cant part of the eigenvalue equation (|) is associated only 
with a small submatrix. 

With this observation, we turn now to calculate the 
main matrix elements of Eq. (||). First, it is straightfor- 
ward to see that 



(1|1) = 1, and 
(1|-1) = (1|±2) = <1|3) = (2|2) = 0. 

Other matrix elements are given by 

(l|-3) = (2|-2)=T 4 = -f - l £ + 0(k 4 



(2|-3) 

(2|3) = |T B | 2 = 2n „ 
where 



k 11 

3072 

- - jl (k s 



rpl 



0(k 5 ), 



(12) 



(13) 



T v = / dxe- 2m[vx+f(x)] . 
Jo 

The above results demonstrate the strong dependence of 
the matrix elements (m\n) on \n\ + \m\, when k < 1. 



k 


Zl 


Z2 


Z3,4 


•25,6 


0.1 


0.035 


-0.037 


-0.003±0.050i 


0.002 ± 0.034i 




0.034 


-0.036 


-0.001 ± 0.035i 




0.2 


0.070 


-0.080 


-0.009±0.098i 


0.011 ± 0.064i 




0.066 


-0.076 


-0.005±0.071z 




0.3 


0.109 


-0.134 


-0.015±0.144i 


0.030 ± 0.088i 




0.094 


-0.12 


-0.01±0.11i 




0.4 


0.153 


-0.196 


-0.018±0.188z 


0.066 ±0.113i 




0.12 


-0.16 


-0.02±0.14i 




0.5 


0.208 


-0.266 


-0.023±0.229i 


0.111 ±0.152i 




0.14 


-0.21 


-0.03±0.18i 





TABLE I. The leading Ruelle resonances of the perturbed 
cat map at various values of the perturbation parameter, k. In 
each case, the top line is the exact numerical value, while the 
bottom line is the result of the variational calculation (ffil). 
The eigenvalue, Zo = 1, associated with the invariant density 
is omitted. 
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FIG. 2. The configuration, in complex plain, of the Ruelle 
resonances of two generic chaotic maps. The right panel cor- 
responds the perturbed cat map at k=0.5, and the left panel 
to the standard map at K = 13. In both cases the coarse 
graining parameter is s — 0.005. 

If we truncate the matrices in @ to 4 x 4 matrices, 
the leading Ruelle resonances are the zeros of the char- 
acteristic polynomial: 



/ 



Det 



o 

T 4 



V Ti - zT 4 





—z 




T, 



4 



—z 




(T 3 2 )* - zTt 
Tt 




These zeros can be calculated exactly, but, having only 
truncated matrices, it is sensible to evaluate them only 
to the leading orders in k. The results are: 



Zl.2 ^ ± 



V8 



z 3A ~ ±i 



k k 2 



(14) 



In Table I we present the numerical (top line) and the 
variational results (bottom line) for the Ruelle resonances 
of the perturbed cat map at various values of the pertur- 
bation parameter k. The numerical calculation is per- 
formed by projecting the equations (|l|) onto a grid of 
90 x 90 sites, and calculating the eigenvalues at s = 
by extrapolation from the interval 0.001 < s < 0.005. It 
has been checked that the results are independent of the 
grid, provided the lattice constant is sufficiently small. 
Comparing the results of Eq. ([u]) to the exact numerical 
values, shows that the variational approach provides the 
order of magnitude and the correct configuration of the 
leading Ruelle resonances. 

To check the generality of this variational approach we 
turn now to consider the standard mapQ as our second 
example. The standard map is defined by 



Vn 



hi = x n + y n 

-l = Vn + g(x n +i) 



(Mod 1) 



(15) 



where 



K 
2^ 



sin(27rcc), 



and K is the stochasticity parameter. The map is inte- 
grable when K = 0, and it becomes increasingly chaotic 
as K is turned on to a large value. The route into chaos 



follows the Kolmogorov-Arnold-Moser scenario of break- 
ing of resonant tori. Therefore, the standard map repre- 
sents a generic mixed system. Yet, for large K, the is- 
lands of stability are tiny and have negligible influence on 
correlation functions of sufficiently smooth observables. 

Applying the same procedure described above, we ob- 
tain that, as before, Eqs. ( |l2| ) are satisfied, but the other 
relevant matrix elements are ( 1 1 — 3) = (2| — 2) = J 2 (K), 
(2 1 -3) = Jf(K), and (2|3) = Ji(K), where J V (K) is the 
Bessel function of integer order. Inserting these results 
into the truncated eigenvalue equation and calculating 
its zeros we obtain: 



21,2 



f 1+ {K)±^r,l(K) + J 2 {K), 
z 3A = V _(K)±^ v 2 _(K) 



where 



V±(K) 



MK), 



J 2 {K)±J 2 (K) 
2J 2 (K) 



(16) 



Thus, to the leading .-approximation in 1/K, the first 
Ruelle resonances areO ±\J 2 (K)\ 1 / 2 and ±i\J 2 (K)\ 1 / 2 . 
As evident from Eqs. (|l6|), z\ <2 diverge whenever K is 
a zero of J 2 (K). At these points our variational ap- 
proach breaks down, but away from them, the results 
are in agreement with the numerical diagonalization, as 
demonstrated in Table II. 

The intriguing feature of the above results is that the 
configuration of the leading Ruelle resonances in both ex- 
amples is similar. The four subleading resonances are lo- 
cated, approximately, at the roots of the equation z 4 = 7 
where 7 characterizes the stochasticity of the map. The 
more stochastic is the map, the smaller is 7. In par- 
ticular, for the perturbed cat map 7 = fc 4 /64 while for 
the standard map 7 = J 2 (K). It is suggestive that this 
behavior is generic to a wide class of chaotic maps. 

In Fig. 2 we depict the positions, in the complex plain, 
of the leading 24 resonances of the perturbed cat map and 
the standard map with a finite coarse graining, s = 0.005. 
This figure suggests that the similarity between the clas- 
sical spectral properties of the map may extend beyond 
the four subleading resonances. However, at this stage 
we do not know how to quantify this similarity. 



K 


21 


22 


23,4 


10 


0.577 


-0.526 


-0.064±0.521i 




0.515 


-0.494 


-0.003±0.505i 


13 


0.617 


-0.561 


-0.002 ± 0.469i 




0.455 


-0.478 


-0.011 ± 0.466i 



TABLE II. The numerical and the variational results for 
the Ruelle resonances of the standard map at K = 10 and 
K — 13. In each case, the top line is the exact numerical 
value, while the bottom line is the result of the variational 
calculation (Eqs. 16). 
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It is natural to ask what are the implications of the 
above results for the spectral statistics of the correspond- 
ing quantum maps? To give a partial answer to this ques- 
tion, we consider the form factor which is the Fourier 
transform of the spectral two-point correlation function. 
Assuming the map to belong to the orthogonal ensemble, 
the semiclassical approximation to the form factor iaJ 



S{n)~2nJ2 



(17) 



a=0 



where n denotes an integer time.— assumed to be much 
smaller than the Heisenberg timell3, and z a are the Ru- 
elle resonances of the corresponding classical map. The 
leading resonance, zq = 1, is associated with the univer- 
sal result of RMT, while higher resonances give nonuni- 
versal contributions. As the map becomes more chaotic 
all \z a \ approach zero except for zq, therefore the non- 
universal corrections to RMT become small. Below we 
show that the configuration of the resonances also plays 
an important role in suppressing the magnitude of the 
nonuniversal contributions. 

Let us assume that z a w for a > 4, and approximate 
the subleading resonances by ±Z and ±iZ, where Z is 
a real positive number smaller than unity. Substituting 
these z a -s in (|l7|) one obtains: 



S(n) fa 2n + 



8nZ n if j is an integer 
otherwise 



From this formula it follows that the nonuniversal correc- 
tions to RMT appear only in powers of Z 4 , rather than 
Z. This is a result of cancellations among the contribu- 
tions of the Ruelle resonances. Thus the magnitude of 
the nonuniversal contribution to the form factor is dic- 
tated both by the absolute values of z a -s, as well as by 
their configuration in the complex plain. It is plausi- 
ble that other chaotic systems exist where the configura- 
tion of resonances lead to an even stronger suppression 
of nonuniversal contribution, e.g. if the subleading res- 
onances are approximately the roots of z = 7, where 
v > 3. 

To summarize, in this paper we have studied the lead- 
ing Ruelle resonances of two maps representing typical 
chaotic behaviors: The perturbed cat map which exhibits 
hard chaos, and the standard map which is a mixed sys- 
tem. Our analytical and numerical results show that, in 
both cases, the configuration of the leading Ruelle, in the 
complex plain, is similar. Numerical studies (e.g. Fig. 2) 
suggest that the similarities in the classical spectrum of 
chaotic maps go beyond the properties of the first four 
subleading resonances. A comprehensive understanding 
of the classical spectral properties of chaotic systems will 
open the possibility for understanding the behavior of 
their quantum counterparts. In particular the weak lo- 
calization mechanisms associated with quantum interfer- 
ence. In this work we show that the configuration of the 
Ruelle resonances may result in a suppression of nonuni- 
versal contribution to the form factor. This mechanism 



of suppression is different from that of diffusive systems 
in which only the magnitude of the Ruelle resonances is 
important. 
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